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Abstract 

We show how to incorporate the lepton sector in a supersymmetric 
theory of flavor based on the discrete flavor group (83)^. Assuming 
that ah possible nonrenormahzable operators are generated at the 
Planck scale, we show that the transformation properties of the leptons 
and of the flavor-symmetry breaking fields are uniquely determined. 
We then demonstrate that the model has a viable phenomenology and 
makes one very striking prediction: the nucleon decays predominantly 
to Kl where I is a first generation lepton. We show that the modes 
n K^^e^ P — ^ K~^iye, and p — > K^e~^ occur at comparable rates, 
and could well be discovered simultaneously at the SuperKamiokande 
experiment. 
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1 Introduction 



The origin of flavor has been a significant puzzle in particle physics since 
the discovery of the muon. The replication of fermion generations and the 
strongly hierarchical pattern of their masses and mixing angles is left un- 
explained in the Standard Model. In most theories that attempt to eluci- 
date the puzzling features of flavor, new symmetries are introduced at mass 
scales that are large compared to the electroweak scale. The new scales may 
be associated with the breaking of flavor symmetries or of a grand unified 
gauge group. In either case, the introduction of a very high scale in the 
theory induces a large radiative correction to the Higgs mass squared, which 
destabilizes the hierarchy between the high scale and the electroweak scale. 
Supersymmetry is the most promising mechanism for avoiding this problem. 
Therefore, it is natural to consider the physics of flavor in the framework of 
the super symmetric standard model (SSM). 

Supersymmetry, however, complicates the problem of flavor by introduc- 
ing a new sector of particles whose masses and mixing angles must also be 
understood. While no superpartner has yet been observed, the acceptable 
spectrum is constrained by low-energy processes. Most notably, a high degree 
of degeneracy is required among the light generation squarks to suppress dan- 
gerous flavor- changing effects unless there is a strong alignment of quark 
and squark eigenstates p|, |^. The challenge in a supersymmetric theory of 
flavor is to simultaneously explain both the suppression of flavor changing ef- 
fects from the scalar sector and the hierarchical pattern of the quark Yukawa 
couplings. Flavor symmetries, spontaneously broken by a hierarchy of vac- 
uum expectation values (vevs), provide an interesting tool for this purpose 
[|, 1^. However, there is considerable freedom in the choice of flavor group 
and symmetry breaking pattern. 

In the recent literature, a number of authors have tried to meet this 
challenge by constructing models of flavor based on Abelian horizontal sym- 
metries, often motivated by superstring theory |Q. However, many of these 
models have problems with large flavor- changing effects The only sur- 
viving explicit Abelian models of which we are aware are those of refer- 
ence 0, which, however, rely on somewhat ad hoc sets of charge assignments 
to achieve alignment. Several authors have considered non- Abelian flavor 
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groups, leading to near degeneracy of the lightest two generation scalars, in- 
cluding SU{2) 0(2) X f/(l) and U{2) §, although the latter appears 
not to solve the flavor changing problem associated with the e parameter of 
K physics. A model has also been proposed in which the three generations 
of fermions are unified into an irreducible multiplet of a non-Abelian discrete 
symmetry This leads to sufficient squark degeneracy to solve the flavor 
changing problem, but a heavy top quark is not guaranteed, and results only 
by assuming that one of three Higgs fields remains light after flavor symmetry 
breaking. 

Two of the authors (LJH and HM) have advocated the use of discrete, 
gauged non-Abelian family symmetries to obtain the desired degree of squark 
degeneracy [l^ : global continuous symmetries are broken by quantum grav- 



itational effects [|TT|, while gauged continuous symmetries may generate D- 
term contributions to the squark and slepton masses that are nonuniversal 
[|T^ , as in the model of reference |Q. It was demonstrated in Ref. that the 
non-Abelian discrete group, (5*3)^, is a promising choice for the flavor symme- 
try group of the SSM. The group S3 has both a doublet (2) and a non-trivial 
singlet representation (1^) into which the three generations of fermions can 
be embedded. In order to construct a viable model, three separate 5*3 fac- 
tors are required, for the left-handed doublet fields Q, and the right-handed 
singlet fields If^ and D'^. The first and second generation fields transform as 
doublets, which ensures the degeneracy among the light generation squarks. 
The third generation fields must then transform as l^s so that the theory 
is free of discrete gauge anomahes. While the group S3 acts identically on 
three objects, the representation structure distinguishes between the gener- 
ations. Thus, it is possible to choose the quantum numbers for the Higgs 
fields so that only the top Yukawa coupling is allowed in the symmetry limit. 
The hierarchical structure of Yukawa matrices can then be understood as a 
consequence of the sequential breaking of the flavor symmetry group. The 
model is appealing on more general grounds since discrete gauge symmetries 
arise naturally in superstring compactifications. 

The model proposed in Ref. [|1^, however, did not address the problem 
of flavor in the lepton sector. The lepton Yukawa matrix is clearly hierar- 
chical, and its eigenvalues are similar in size to those of the down quarks. 
In addition, a high degree of degeneracy between first and second generation 
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sleptons (or an alignment between slepton and lepton mass eigenstates) is 
required to suppress dangerous lepton flavor-violating processes. As far as 
this point is concerned, it is reasonable to expect that the flavor structure 
proposed in Ref. for the quarks should work equally well when applied to 
the lepton sector. However, there is no guarantee that the flavor symmetry 
will provide an adequate suppression of the operators which mediate pro- 
ton decay. Recall that there are two possible sources of proton decay in the 
SSM: i?-parity violating dimension-four operators (like QDL or UDD), and 
nonrenormalizable, dimension-five operators (like QQQL or UUDE) which 
are likely to be generated at the Planck scale. In the first case, the coeffi- 
cients of the i?-parity violating operators are forced to be extremely small 
by the nonobservation of proton decay. The operators ciQsL and C2uds, for 
example, are constrained such that C1C2 < 0(10^^^). In this paper we sim- 
ply assume that i?-parity is an exact symmetry, and these operators are not 
present.]] Assuming that quantum gravity effects violate global symmetries, 
however, non-renormalizable operators that conserve i?-parity, but violate 
baryon and lepton number, are presumably generated at the Planck scale. 
The coefficients of the operators {QiQi,2)Q2Li/ Mp are constrained to be no 
larger than 0(A®) with A ^ 0.22. In the absence of other mechanisms to 
eliminate these operators,! an adequate theory of flavor must explain why 
they are sufficiently suppressed after flavor symmetry-breaking effects are 
taken into account . 

In this article, we extend the {S^Y niodel to the lepton sector. First, 
we require that the leptons transform under the same 5*3 flavor groups as 
the quarks. This is the simplest choice given that the ordinary Higgs fields 
transform nontrivially under and 5*1^. We select the transformation prop- 
erties of the lepton fields under (6*3)^ so that we obtain the greatest similarity 
between the lepton and down-quark Yukawa matrices. The assignment must 
also forbid all dangerous dimension-five operators in the (5*3)^ symmetry 

^Since matter parity is non- anomalous with the minimal SSM (MSSM) particle content, 
it can be considered as a discrete gauge symmetry, and is hence preserved by quantum 
gravity effects. i?-parity is a product of matter parity and a 27r rotation in the local 
Lorentz frame. 

^Yov example, we could impose a Peccei-Quinn symmetry [jl^ , discrete symmetries 
ill, or gauge U(1)b @. 
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limit. We will see that the transformation properties of the lepton fields are 
uniquely fixed by these requirements. We then argue that the fundamental 
sources of fiavor symmetry breaking are gauge singlet fields that transform 
in the same way as the irreducible "blocks" of the quark Yukawa matrices. 
We will call these fields 'fiavons' below. We show from a general operator 
analysis that models involving fiavons with simpler transformation proper- 
ties can all be excluded, if flavor physics originates at the Planck scale. With 
the flavor symmetry breaking originating only from the Yukawa matrices, we 
consider the contributions to lepton flavor violation and proton decay. We 
show that the model is consistent with the current experimental bounds. In 
addition, we show that the dominant proton decay modes in the (Ss)^ model 
are of the form p ^ Kl, where I is a first generation lepton. This is never the 
case in either supersymmetric or non-supersymmetric grand unifled theories. 
The prediction of these rather unique modes is exciting since the total decay 
rate is likely to be within the reach of the SuperKamiokande experiment. 

In the next section we make the assumptions of our analysis explicit. In 
section 3, we introduce an economical form for the flavons (p, which provide 
an adequate description of both lepton and quark mass matrices via dimen- 
sion flve interactions. We use flavor-changing and baryon number violating 
phenomenology to demonstrate that this choice is the simplest possible ac- 
ceptable form for the fiavons in section 5. This phenomenology is studied in 
much further detail in section 6, and conclusions are drawn in section 7. 

2 The Framework 

In this paper we construct a description of flavor, for both quarks and leptons, 

based on a fiavor group (<S'3)^, spontaneously broken by a set of fiavon fields, 
(f). We write down an effective theory beneath the Planck scale in which the 
gauge symmetry is Gsm = SU (3) x SU (2) x U (1) and R parity is imposed. 
The field content is that of the minimal supersymmetric model, together with 
a set of fiavon fields, which are necessarily all neutral under the gauge group 
since they take vevs much larger than the weak scale. The theory contains 
the most general F and D terms consistent with Gsm x [Ss)^ x Rp with all 
interactions scaled by the appropriate powers of Mpi and all dimensionless 
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coefficients of order unity. The single exception to this is the absence of a 
Planck scale mass for the Higgs doublets and for the flavons. In addition the 
theory is taken to possess supersymmetry breaking interactions which are 
the most general soft operators consistent with the symmetries of the theory. 
In particular, no universality assumptions are made to relate otherwise free 
parameters. At the renormalizable level, the theory is remarkably simple: 
the only F terms are the Yukawa coupling for the top quark, and possi- 
ble trilinears amongst the flavon fields. The only supersymmetry breaking 
terms are: the three gaugino masses, a trilinear scalar interaction involving 
the top squarks coupled to a Higgs doublet, scalar masses for the flavons, 
the Higgs, and for the squarks and sleptons, which are diagonal in flavor 
space with the structure (m^,m^,m|). Supersymmetric non-renormalizable 
F terms lead to the Yukawa matrices becoming functions of (p/Mpi, while 
the non-renormalizable F and D terms, which contain supersymmetry break- 
ing spurious, similarly lead to the trilinear A terms and the scalar masses 
becoming functions of (p/Mpi. 



3 The Basic (2,2) Model 



In the (Ss) model of Ref. [|T0|, the quark chiral superfields Q, U, and D are 



assigned to 2 + 1a representations of S^, and S^, respectively. The ffist 
two generation fields are embedded in the doublet, for the reasons described 
in the Introduction. The Higgs fields both transform as (1^, 1a, l5)'s, so 
that the top quark Yukawa coupling is invariant under the flavor symmetry 
group. The transformation properties of the Yukawa matrices are: 



Y, 



u 



(2, 


2, 


I5) 


(2,15,15) 


(l5 


,2 


I5) 


(15,15,15) 



Y, 



D 



(2,1a, 2) 


(2, 1a, 1a) 


(l5,lA,2) 


(15,1a,1a) 



(1) 

where we use the notation 2 = 2 (g) 1^,|. Note that these matrices involve at 
most 7 irreducible multiplets of (5*3)'^. In Ref. [0, (S^)^ was broken by only 
four types of flavons: 0(2, 15, 15), 0(2, 2, 15), 0(1^, 1a, 1a), and 0(2, 1^, 2), 
the minimal number which allows us to obtain realistic masses and mixings 



§2 = (a, 6) is equivalent to 2 = (6, —a). 
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breaking 


^3 


CU cD 
^3 "-^3 


htVcb 




cU cD 

3 3 


hb 




(^3^ X S^)/Z, 


htVub 


nothing 


(s!^ X s^)/Z2 


he, hcX 


nothing 


{Z^ X S^)/Z, 


hs, hgX 


nothing 


ryU,D 

^2 


hu, hd 


nothing 


nothing nothing 



Table 1: The sequential symmetry breaking pattern of (5*3)^ symmetry. The 
degeneracy between d and s is kept until the non-Abelian group (Z^f x S^)/Z2 
is broken. Z^'^ symmetry remains unbroken until the last stage, which keeps 
up- and down- quarks massless. 



ni: 





' hu 


hcX 


-htVub 


\ 




( hd 


hsX 


\ 


Yu = 





he 


-htVeb 




, Yd = 





hs 







I 





ht 


) 




I 





hb ) 



with A ~ 0.22. We will refer to this scenario as the basic (2,2) model below, 
because it involves flavons which transform as doublets under two of the 
5*3 factors simultaneously. Note that it is not absolutely necessary that we 
generate Veb and Vub by rotations in the up sector, as we have indicated 
above. We can instead generate Veb and Vub in the down sector (by assigning 
appropriate values to the (1,3) and (2,3) elements in Yd) using the breaking 
parameter 0(2, 1a, 1a) rather than 0(2, Is, Is)- This choice generally gives 
us much weaker phenomenological constraints, as we will see later. In section 
5 we study whether one can construct 0(2, 2, Ig) and 0(2,1^,2) breaking 
parameters from products of flavons that transform as doublets under only 
one 5*3, and find it is not possible within the framework specified in the 
previous section. 

The form of the Yukawa matrices presented above can be understood as 
a consequence of the sequential breaking of the flavor symmetry, as shown in 
Table |I|. Note that it is necessary to have different flavons associated with 
each step in the sequence, since all the components in a single irreducible 
multiplet typically become heavy at the same time unless a fine-tuning is 
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done. For instance, we assume that one chiral superfield transforming as 
0(2, 15,15) acquires a vacuum expectation value (VEV) in its Vi component, 
which breaks to Z2 and generates htVcb- A different chiral superfield 
transforming in the same way, remains light at this stage, but acquires a VEV 
in its V2 component at lower scale to break Z2 and generate htVub- Since 
is completely broken at this stage, a single chiral superfield 0(2, 1^, 2) is split 
into two (1^, 2)s under the remaining [S^ x S^)/Z2 symmetry; both acquire 
VEVs in their Vi components to break this symmetry down to Since 
hs and hsVcd are generated at the same stage of the symmetry breaking, there 
is a natural reason why the Cabibbo angle is rather large. The final stage 
of breaking is done by another 0(2, 1a, 2) to generate ha (and by another 
0(2, 2, 15) to generate hu in the up sector). Therefore, the 2x2 block in Yd 
has the structure Yd = + ^2, where 

Note that Yi preserves Z^'^ symmetry. Of course the other elements in Y2 
can be non-vanishing, but are expected to be of order hd or less, and are 
irrelevant for our purposes. Similarly, the 2x2 block in Yu is given by 
Yu = Y/ + Y2', where 

yi- 1 !! ""'f^^ h yi= I : I- (4) 




h 

Note that a and a' are order one constants, with a — a' — 1. 



4 Incorporating Lepton Sector 

The Higgs fields in the {S^Y model transform nontrivially under the fiavor 
symmetry group. Since the lepton fields acquire their masses in the MSSM 
from the same Higgs fields as the quarks, the leptons should transform under 
the same (5*3)^ fiavor symmetry. We are led by three principles in determining 
the precise transformation properties of the lepton fields: 

1. We do not allow any new fiavor symmetries (e.g. new 6*3 factors) that 
arise only in the lepton sector. The only fiavor symmetry in the theory 

IS iS*2 ^ ^S*^ X ^2 . 
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2. We assign the transformation properties of the lepton fields so that the 
charged lepton Yukawa matrix is similar to that of the down quarks. 
This choice is suggested by the phenomenology. 

3. We require that the most dangerous dimension- five operator that con- 
tributes to proton decay, {QQ){QL), is forbidden in the (S^)^ symmetry 
limit. 

As we will see below, these principles are sufficient to completely determine 
the transformation properties of the lepton fields. 

Let us first consider the consequences of the first two conditions. The 
down-quark Yukawa matrix is a coupling between the left-handed quark 
fields Q ~ (1^ + 2, Is, Is) and the right-handed down quark fields D ~ 
(I5, Is, 1a + 2). We know that the Yukawa matrix of the charged leptons is 
quite similar to that of the down quarks, up to factors of order three [0] at 
high scales: 

nib ~ rur, rus ^ -m^, rrid ~ Sm^. (5) 

Therefore, we look for an assignment of lepton transformation properties 
that leads automatically to this observed similarity. There are only two 
possibilities: 





s? 




sF 








qU 




L 


1a + 2 


Is 


Is 


or 


L 


Is 


Is 


1a + 2 


E 


Is 


Is 


1a + 2 




E 


1a + 2 


Is 


I5 



(6) 



The third condition above allows us to distinguish between these two al- 
ternatives. In the first assignment, the operator {QiQi){QjLj) is allowed 
by the (83)^ symmetry, and we have proton decay at an unacceptable rate. 
Therefore, only the second assignment in Eq. satisfies all three criteria 
listed above. We could have obtained the same conclusion by considering the 
UU DE operators as well. 

The remaining question that we need to answer is how the factors of 
three in enter in the Yukawa matrices. One plausible explanation is that 
they originate from fiuctuations in the order one coefficients that multiply 
the (5*3)^ breaking parameters which generate the quark and lepton Yukawa 
matrices. As discussed in the previous section, it is quite likely that the 
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2x2 block in Yd is generated by two 0(2, 1a, 2) breaking parameters. An 
acceptable lepton Yukawa matrix is obtained by allowing coefficients in both 
breaking parameters to deviate from unity by a factor of three. Throughout 
this paper we take Yj = 3Yi + ^¥2. 

5 Uniqueness 

In the basic (2,2) model introduced in the previous section, the quark Yukawa 
matrices were taken as the only sources of flavor symmetry breaking. Thus, 
the Higgs fields that spontaneously break the flavor symmetry group in the 
corresponding high energy theory would come in exactly four representations 
of (Ssf: (2, 1a, 2), {Is, U, U), (2, 2, and (2, 1^, These correspond 
to the various blocks of the quark Yukawa matrices. We have assumed that 
no other representations are involved in (S^)^ breaking in the full theory. 

The issue that remains to be addressed is whether this picture of the 
high energy theory is overly restrictive. There are 3^ — 1 = 26 nontrivial 
representations of (83)^ that we could have used had we tried to build an 
adequate high-energy theory of flavor symmetry breaking at the start. There 
is no reason to assume a priori that a model cannot be constructed with a 
different, more fundamental set of symmetry breaking parameters. What we 
will show in this section is that there are in fact no simple models of (5*3)^ 
breaking involving symmetry breaking parameters that are more fundamental 
than the ones adopted in the basic (2,2) model. We will systematically 
exclude all the reasonable alternatives. While the basic (2,2) model served 
as an existence proof for a successful {S^Y model in Ref. |jlO|, we will show 
here that the choice of symmetry breaking parameters in this model is in fact 
unique. 

We first would like to consider the class of models in which there are no 
"(2,2)" representations, i.e., there are no fields that transform as doublets 
under more than one S3 group at a time. As a starting point, let us consider 
a simple toy model that illustrates the phenomenological problems common 
to models of this type. The flavon content is 

~ (2, 15, Is) Hu {ls,2,ls) Hn--{ls,ls,2) (7) 
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where i = 1 . . .2 labels two distinct doublets, and 



Xi ~ (I5, 1a, I5) X2 ~ (I5, 1a, 1a) X3 ~ (1a, 1a, I5) X4 ~ (1a, 1a, 1a) 

(8) 

The doublet H fields were chosen to have the simplest quantum number 
assignments possible. The quantum number assignments of the x fields were 
chosen for a variety of reasons, that will become clear in context below. 
When the fiavon fields acquire vevs, the various blocks of the quark mass 
matrices are generated from higher dimension operators. The way in which 
this model fails is instructive and will simplify the discussion of the other 
models to follow, so we will proceed in some detail. 

The two-by-two quark Yukawa matrices in this model are obtained by 
taking the products (T2 i^gi^^crj and a2HQHDXu-i for the up and down sectors 
respectively. Given this construction, we require two Hq fields in order to 
assure a nonvanishing Cabibbo angle. Let us denote the ratio of the vevs of 
the H and x fields to an appropriate cutoff scale by e and 5, respectively. If 
we take the Hq and Hd to be of the form 



1 

A 



Hd — 



A 
1 



(9) 



then we obtain the two-by-two down Yukawa matrix 



A2 A 
A 1 



(10) 



where A ~ 0.22 is the Cabibbo angle. By setting the combination tQe^Si ~ 
A^, we obtain the correct strange quark Yukawa coupling, assuming tan/3 ~ 
1. 

The down quark Yukawa coupling has not yet been generated, however, 
because the matrix above has a zero eigenvalue. This is a generic feature 
of all models in which the down-strange Yukawa matrix is formed by taking 
the product of two doublets. Of course, the same problem arises in the up- 
charm Yukawa matrix as well. The way in which the up and down Yukawa 
couplings are generated is through other operators that arc of higher order 
in the symmetry breaking. Notice that we can obtain a correction to the 
down-strange Yukawa matrix via the operator X4ifQif|^aJ which is of the 
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form 



1 A 
A A2 



(11) 



If we take 54/^2 ~ A^, then we obtain a correction to the (1,1) component of 
(p^Ol). This is sufficient to hft the zero eigenvalue at order eQe£)(54 ~ A^, which 
is of the correct order in A to generate the down Yukawa couphng. A similar 
mechanism occurs in the up sector as well. If we take 



(0 



eg 



1 

A 



(12) 



then the up-charm Yukawa matrix a2HQH^aJ is given by 



A2 A 
A 1 



We take egeu ~ A^ in order to generate the correct charm quark Yukawa cou- 
pling. The zero eigenvalue in this matrix is lifted via the operator xsHqH^, 
which is of the form 

' 1 A 



A A^ 



If we take ^3 ~ A^, then this corrects the (1,1) entry of the ffist operator, 
lifting the zero eigenvalue at order A^, as desired. 

The inescapable problems with the model above arise from flavor changing 
neutral currents and proton decay considerations. In the ffist case, the off- 
diagonal elements in the squark mass matrices can be constructed at ffist 
order in the symmetry breaking e parameters via the '2^ invariant' 



iQ*' asQ)H, - iQ*^aiQ)H2 



(13) 



Here Q is the squark doublet that transforms under the same 5*3 group as H, 
and we use subscripts to denote the components of the H doublet. The result- 
ing off-diagonal elements are constrained by flavor changing neutral currents 



processes, such that {Hq2) < 0.05, (-ff_D2) < 0.05, and {Hq2){H£)2) < 
0.006. Given the form of the H fields described above, this implies that 



egeijA < 0.006. However, we saw earlier that the strange quark Yukawa 
coupling hg ~ cqCd^i is of order A^. Hence we require 61 ~ 3, which con- 
tradicts our assumption that all the e's and 5's are small symmetry breaking 
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parameters. In other words, for 6i < 1, we cannot generate a large enough 
strange quark Yukawa couphng if we are to simultaneously satisfy the FCNC 
constraints. 

The second problem is that there are operators in this model that con- 
tribute to proton decay at an unacceptable level. We can construct the rep- 
resentation (2,15,2) at order egtr), which contributes to proton decay via 
the operator -^{QQQL) at the same order. This forces us to take ege^ < A^, 
which again makes it impossible to generate a large enough strange quark 
Yukawa coupling. 

What we have found is that it is impossible to generate large enough 
down quark Yukawa couplings when the H fields have the simple transfor- 
mation properties described above. One way to remedy this problem is to 
allow these fields to transform under more than one Sz group, so that the 
contributions to the off-diagonal elements of the corresponding squark mass 
matrices occur at higher order, and dangerous proton decay operators are 
sufficiently suppressed. The minimal modification of the model above with 
this property is a model in which some or all of the H fields transform as IaS 
under an additional (or both remaining) 5*3 groups. Since it is possible in 
this case to build the down-strange Yukawa matrix from the product of Hq 
and Hd alone, we will restrict our discussion to models in which the both 
the down-strange and up-charm Yukawa matrices are generated at second 
order in the symmetry-breaking e parameters. Any model in which these 
matrices are generated at higher order in the symmetry breaking will involve 
the representations described below as composite operators, and at the very 
least will be subject to the same constraints. There are a finite number of 
possibilities for the type of model of interest, and we will now outline why 
each fails: 

case 1: Hd ~ (lyi. Is, 2). This representation contributes to proton decay 
through the operator j^{QQQL) at order X^ec,. Assuming that none of the 
e's arc larger than order A, then eo > A^ so that we can generate a large 
enough strange quark Yukawa coupling. Hence, the coefficient of the proton 
decay operator is order A^ or greater. This leads to an enhancement in the 
proton decay rate by four orders of magnitude over that of the basic (2,2) 
model, and this possibility is excluded. 

case 2: Hd ~ (Is, 1^4, 2). To construct the upper-left two by two block 
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of the down quark mass matrix using this representation we must also have 
Hq^ ~ (2,15,15). (The alternative choice Hg^ ~ (2, 15,1^), for any i, is 
excluded by the proton decay constraints, as we will discuss later.) Notice 
that we generate the (1,3) and (2,3) entries of the down quark matrix from 
the product of Hq and hb. Thus we require eg to be order or smaller if 
Vub and Vcb are not to be unacceptably large. This forces us to take of 
order or larger if we are to generate an adequate strange quark Yukawa 
coupling. Now the problem arises because Hd also contributes to the (3,1) 
and (3,2) entries of the down quark Yukawa matrix at order eo, which we 
have just argued is of order A^ or larger. This gives us at least an order A 
rotation on the right-handed down quark fields between the first and third 
generations. In the basis where the quark Yukawa matrices are diagonal, this 
yields a (3,1) entry in the right-handed down squark mass matrix of order 
AMi (Recall that the diagonal components of the squark mass matrices Mi 
and M3 are unconstrained.) This is in marginal disagreement with the the 
bound {m^)f^/Mf < 0.1 from the B^-^ constraint. A more serious problem 
arises because there is also an order 1 rotation between the right-handed b and 
s fields, which spoils the degeneracy between the second and first generation 
squarks; since there is an order A rotation on the right-handed squarks of 
the first two generations in this model, we obtain a (5^^) 12 of order A, which 
exceeds the bound (5^^)i2 < 0.05. 

case 3: Hn ~ (1^, 1a, 2) In models that include this representation we 
must also have Hq^ ~ (2, 15, 15). (Again, the alternative Hq^ ~ (2, 15, 1^), 
for any i, is excluded by proton decay constraints, as discussed below.) Hq 
contributes to the (1,3) and (2,3) entries of both the up and down Yukawa 
matrices at order eg and tghb respectively, but always with the largest com- 
ponent in the (3,1) entry. Thus, this model predicts Vub > Vcb, and is there- 
fore excluded. 

Now that we have established that H^, ~ (15,15,2) is the only viable 
representation for Hd (that has no more than 2 degrees of freedom) we are 
forced to modify the Hq if we are going to evade the problems of the simple 
model presented at the beginning of this section. 

case 4-' Hq ~ (2, 1a, 1a)- This representation contributes at order eg to 
the (1,3) and (2,3) entries of the down quark Yukawa matrix. Thus, eg can 
be no larger than order A^ if Vcb and Vub are to be small enough. In this case. 
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we cannot generate a large enough strange quark Yukawa coupling. 

case 5: Hq ~ {2, Is, 1a)- This contributes to proton decay via the 
■^{QQQL) operator at order eg and is immediately excluded. 

case 6: Hq ~ (2,1^,15). We can construct the spurion (2, 15, 1^) at 
order hi^eq ~ A^e^, which contributes to proton decay via -^{QQQL) at the 
same order. This forces us to take eg ~ A'^ or smaller, which in turn tells 
us that cd ~ A, in order to generate a large enough strange quark Yukawa 
coupling. However, this is excluded by the flavor changing neutral current 
constraints for Ho ~ (Is, Is, 2). 

The discussion above forces us to consider representations for the if's 
that have four degrees of freedom. Let us consider the following possibilities: 

case 7: Hd ~ (2, Is, 2) . We could imagine constructing a model in 
which the down two by two yukawa matrix is a product of (2, Is, 2) and h},. 
However, Hd now gives us a contribution to proton decay at order e^, and 
hence this alternative is excluded. 

case 8: Other Models without a (2, 1^,2). We still might hope to con- 
struct models without a (2, 1^2) if we can generate the upper two by two 
block of the down quark Yukawa matrix at second order in the symmetry 
breaking. In such models, we require that the upper two by two block of 
the up quark Yukawa matrix be generated at second order as well. This 
restricts us to models with the representations (a) (1^, 2, 2) and (2, 2, 1^) 
or (b) (Is, 2,2) and (2,2,1^). In both cases, the down two by two block 
is generated by taking the product of the two spurious, while the up matrix 
is generated by taking the product of (2, 2, 1a) and h}j. Of course, we must 
also introduce a (2, Is, Is) if we are to generate sufficient Vuh and Vch- The 
problem with these models is that we can also construct a (2, Is, 2) by tak- 
ing a product of the same spurious that generate the strange quark Yukawa 
coupling. Hence, there is a contribution to proton decay at order A^, and 
these models are excluded. 

We have seen above that we must have a (2, 1^, 2) as a fundamental field 
in the model, and that the down and strange Yukawa couplings are gener- 
ated at first order in the symmetry breaking. If we require that the up quark 
two by two block also be generated at first order (so that the symmetry 
breaking is of a comparable magnitude) then we must introduce a funda- 
mental (2, 2, 15). Since we cannot generate large enough values for Vuh and 
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Vcb by taking products of these representations alone, we must also introduce 
a (2, Is, Is). Finally, we generate hh most economically with a (l^, 1a, 1a)- 
Thus, we have arrived uniquely at the basic (2,2) model described in sections 
3 and 4. 



6 Phenomenology 

In the previous sections, we showed that the transformation properties of 
both the lepton and flavon fields are uniquely determined in the {S^Y model, 
if all possible nonrenormalizable operators are generated at the Planck scale. 
With these results at hand, we may now consider lepton flavor violation and 
proton decay, two topics that depend crucially on the extension of the model 
to the lepton sector. We will also consider the bounds from CP-violating 
processes, which were not covered in the original paper. We demonstrate 
that the model is indeed viable phenomenologically. (Whenever we discuss 
numerical estimates, we take a representative choice tan/? ^ 2.) In addition, 
we show that the model predicts the dominance of the Ke proton decay 
mode over the Kfi or vre mode, which is in sharp contrast to the situation in 
supersymmetric or non-supersymmetric GUTs. 



6.1 Lepton Flavor Violation 

The strongest constraint on lepton flavor violation comes from the non- 
observation of the fi —>■ ej decay mode. In our model, the contribution 
of the off-diagonal term in the purely left-handed slepton mass matrix (the 
LL matrix) is small enough (~ h^X ~ 1 x 10~^) to avoid the experimental 
constraint for any value of mj above the LEP bound. The stringent limits 
come from the purely right-handed slepton mass matrix (RR) and the left- 
right (LR) matrix, which we discuss in this section. We follow the notation 
of Barbieri, Hall and Strumia |T3 below. For simplicity, we work in the 



approximation where the exchanged neutralino is a pure bino state. 

The one-loop slepton and bino exchange diagram that picks up the off- 
diagonal (2,1) component in LR mass matrix generates the operator 

O = ^F^iml, ml, Ml) eRta^yLF,. , (14) 
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where is the electromagnetic field strength and 

G2iml,Mll-G2{mlM^ 

with 

M 



2(m2 - M2)3 
The decay width is given by 



- - 2m^M^ In 



M2 



(16) 



Tif, ^ e^) = -mllF^W (17) 

and the bound Br(/i ^ 67) < 4.9 x 10^^^ implies IF2I < 2.6 x 10"^^ GeV^^ 
In order to compare this bound to the prediction of our model, let us take 
mji = rriL = m = 300 GeV and Mi = 100 GeV as a representative case. 
Note that F2 is at its maximum around Mi ~ m/2, and our choice gives 
almost the largest possible F2 for fixed m. We obtain 

< 1.0 X 10-^ (18) 

for this choice of parameters. In our model, the (2,2) and (1,2) elements 
in Yi belong to the same irreducible multiplet, and diagonalization of Yi 
also diagonalizes LR mass matrix at 0{hsX). The term which may not be 
simultaneously diagonalizable comes from the piece Y2 ~ h^, and hence 

(^D2i ~ m^XA, (19) 

where is the down quark mass evaluated at the Planck scale ~ 
10 MeV/3, and A is a typical trilinear coupling. If we take A ~ 100 GeV, 
then {mj^p()2i/m^ ~ 0.8 x 10~^ and the constraint ( [T^ ) is easily satisfied. 
Since we have consistently allowed an order 3 ambiguity in estimating the 
coefficients of various operators, and we do not know the precise value of m^, 
the actual constraint is weaker than the one considered above. The (1,2) 
element in the LR mass matrix contributes in exactly the same way as the 



(2,1) entry, except that the chiralities of the electron and muon in eq. (|lj) 
are fiipped. Hence, the (1,2) element is subject to the same constraint, which 
again is clearly satisfied in our model. 
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The RR mass matrix also contributes to the operator in (p!^. In this 
case, the function F2 is given by 

= ^^.(^L)i2 ''g^, '\ (20) 

where 

dGiim^.Ml) 1 



Sm^ 6 {m? — M^)^ 

(m^ - 9 M2 - 9 m2 + 17 + (18 + 6 M^) log(f^)) . (21) 

This is a monotonically decreasing function of M for fixed m. For the bino 
and slepton masses chosen earher, we obtain the bound 

{mRR)i2 ^ Q_Q23, (22) 



rn? 

while in our model 

(^Di2M' ^ /iiKfeA ~ 0.009 . (23) 

The bound on the (1,2) element of the RR matrix is easily satisfied in our 
model. Note had we chosen the option of generating Vcb and Vuh in the down 
sector, as discussed in section 3, we would have obtained a much smaller off- 
diagonal element {Tn\p^i2- In this case the breaking parameter transforms 
as 0(2, 1^, 1a) rather than 0(2, 15, 15), and does not contribute to the RR 
mass matrix at the first order in 0. The leading contribution is then 

(^L)i2M' - hlV,b^ ~ 3 X 10-*^ , (24) 

and the constraint ( ^2]) is again easily satisfied. 

Finally, there is another contribution to /i — 67 from the mixing to the 
third generation sleptons. Since the third generation scalars can be non- 
degenerate with the first and second generation ones, GIM cancellation does 
not occur. This is similar to the situation in the minimal SO(IO) model 
|l^], where the third generation sleptons are much lighter than the oth- 
ers, while the rotations to the quark mass eigenstate basis are CKM-like. 
Recall that in the minimal SO(IO) model, the flavor changing factors in the 
amplitudes are VtdVts, whereas in our model they are VubVcb, which is about 
a factor of three smaller. 

To summarize, the model is clearly consistent with the experimental 
bounds on lepton flavor violation, and predicts /i — >• 67 at a rate just beyond 
the current limit. 
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6.2 CP Violation 



If we allow arbitrary phases in the symmetry breaking parameters, the squark 
mass matrices can have complex elements which contribute to CP violating 
effects. The most stringent limit on these phases comes from e', and has been 



studied by Gabrielli, Masiero and Silvestrini [^. If we consider the most 
extreme case imaginable, where all the off-diagonal elements of the down 
squark mass matrices in our model are purely imaginary, then we find 



Im 



Im 



Im 



[m 









12 


IV? 










)l2 





/iwA ^ 9 X 10^ 



htVcbX ^ 0.009, 



/i?A ~ 1 X 10 



-7 



(25) 
(26) 
(27) 



in the basis where the quark Yukawa matrices are diagonal. The constraints 
on these elements are 



Im 



m? 



< 2 X 10" 



\ 



Im 



[m 



ll,rr)i'2 
m? 



< 3 X 10" 



\ 



Im 



LL)12 



< 2 X lo- 



rn' 



m' 



(28) 
(29) 
(30) 



for rriq ^ rrig ~ 500 GeV. One can see that all constraints are easily satisfied. 
In addition, there are a number of factors that make the actual bounds 
on our model weaker: (1) This analysis is valid only up to the unknown 
factors of order one that multiply the symmetry breaking operators. (2) The 
renormalization group running of soft SUSY breaking masses always tends to 
make the diagonal elements in the LL, RR, and LR matrices larger at lower 
energies, so that the true constraints on our model are generally weaker than 
those given above. (3) The elements Vch and Vub could instead be generated 
in the down sector, in which case {171^^)12/17? ^ /i^Vc^A, and the SUSY 
contribution to e' becomes negligible. (4) We have no reason to expect that 
all the off-diagonal elements of the squark mass matrices in our model are 
purely imaginary, as we have assumed to obtain these bounds. 
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operator 


{S-^Y representation 


largest coefficient 


{ r\ r\ \ ( r\ t \ 


(2, Is, 1a) 


hshd or hhcK 


[QsQsjiQtLi) 


(2, Is, 2) 




[QiQijiQsLi) 


(2, Is, 2) 


t. t. 


( n n \ (n r \ 


I-La, is, 


h h h \T/ 

ntiibfisAVcb 


iQ^Q^)iQ3L,) 


(2, Is, 1a) 


hshd or hbhchu 


{Q^Q^)iQ3L3) 


(1a, Is, 1a) 


hshd 


{Q^Q^){QiL3) 


(2, Is, 1a) 


hshd or hhchu 


{Q^Q^{Q^L^) 


(2, Is, 2) 


hbhs 



Table 2: A complete list of possible baryon-number violating dimension-five 
operators involving doublet fields. All other operators vanish because of sym- 
metry reasons. Qi and Lj generically refer to either ffist- or second-generation 
fields, while and L3 refer to third-generation fields. For operators that are 
multiplets under the (5*3)^ symmetry, the coefficient of the largest component 
is given. 

As in the case of — >■ ej, there is also a contribution to e' from the third 
generation squarks, similar to the minimal SO(IO) case [^. This falls within 
an acceptable range. Hence, it is reasonable to conclude that the model is 
consistent with the observed CP violating phenomenology. 

6.3 Proton Decay 

Since we have assumed throughout this paper that all possible nonrenor- 
malizable operators are generated at the Planck scale, the task of studying 
proton decay in our model is a simple one. We ffist write down all the possible 
dimension-five operators that contribute to proton decay and identify their 
transformation properties under (S'a)^. The coefficients can be estimated as 
the product of Yukawa couplings that will produce the desired symmetry 
breaking effect. A list of possible operators and their coefficients is given in 
Tables | and |. 

As we can see from the tables, the most important operator involving 
left-handed fields is {QiQi){QiLi) /M^, where = Mpi/y/Wn is the reduced 
Planck mass, and where parentheses indicate a contraction of SU(2) indices. 
This operator transforms as a (2, Is, 2) under (S'3)^, and therefore has a coef- 
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operator 


(os) representation 


largest coerncient 


tUiOT 


(1^,2,1^) 


7 7 7 T 7 

hthbh^Vub 


tuihci 


(2,2,1a) 


7 7 

hhhc 


I 7 

lUidiT 


(1a, 2, 2) 


1- 1- \ 

nciisA 






riciis 


uchr 


(1a, 1a, 1a) 


h^hbVcbVub 


ucbci 


(2, 1a, 1a) 


hthbVcb 


ucdiT 


(1a, 1a, 2) 


hthgVcbX 




(2, 1a, 2) 


hs 



Table 3: A complete list of possible baryon- number violating dimension- 
five operators involving singlet fields. All other operators vanish because of 
symmetry reasons. -Uj, dj, and Cj generically refer to either first- or second- 
generation fields, while t, h and r refer to the third-generation fields. The 
only possible product of the form -Uj-Uj is uc because of the anti-symmetry in 
the color indices. For operators that are multiplets under the {S^)^ symmetry, 
the coefficient of the largest component is given. 

ficient of order hbhg. The two other operators that have the same coefficient, 
{QzQz){QiLi) and {QiQi){QzLi), each involve third generation fields, and 
therefore contribute to proton decay at a rate that is further suppressed by 
small CKM angles. The coefficients for all four components of the leading 
operator are given by 

O = ^^[hs{Q2Qi){QiL^)-hs\{QiQi){Q2L^) 

-0{ha){Q2Q2){QiL2) + ha{QiQ^){Q2L2)\ ■ (31) 

The overall coefficient c is a number of 0{1), and the factor of 1/2 has been 
included to compensate for the symmetry of each term under the interchange 
of either the two QiS or two Q2^s. The striking feature of this multiplet 
is that the operators involving first generation lepton fields Li have larger 
coefficients than those involving second generation fields L2. 

The only operator involving right-handed fields that appears potentially 
dangerous is hs{uc){se — Xsfi). ffowever, the contribution of this operator 
to proton decay is in fact negligible. Since all fields are right-handed, their 
flavor cannot change via l^-exchange. A c squark can change to a m squark 
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Qi, El 


O2, E2 Q:i, Ea 


Ui U2 Ua 


Di, Li Da, L2 Da, La 


Hu Ho 


^2 


+ + + 


+ - - 


+ - - 





Table 4: Charge assignments of MSSM fields under the Z2 ' symmetry. 



via a flavor off-diagonal element of the RR squark mass matrix, and the 
resulting operator can be dressed by a gluino. Since the off-diagonal entry is 
("^iji?)i2/'^^ — h^X, the effective coefficient is smaller than hghc- Thus, this 
operator is negligibly small compared to the leading left-handed operator 
discussed above. The same can be shown for all the remaining operators in 
Table 

What we have concluded based on the leading operator is that our model 
favors proton decay to u^. and e over decay to and fi. This result is in 
striking contrast to the situation in grand unified theories, and provides a 
counterexample to a common theoretical prejudice. Usually it is assumed 
that the operator involving Li should be most strongly suppressed because 
it only violates the electron's U(l) chiral symmetry through the tiny electron 
Yukawa coupling. Nevertheless, the operator involving the electron in our 
model is larger than the operator involving the muon. 

This unusual feature can easily be understood by considering the residual 
^2'^ symmetry that is present when the Yukawa couplings of first generation 
fields are set to zero. Suppose that the all the breaking parameters except 
hu and are present. This leaves a 2'^'^ symmetry, where the fields have 
the charge assignments given in the Table. |^. 2'^^'^ is a diagonal subgroup 
of Z2 and What we learn from this table is that the first generation L 
field is even under this Z2, while the second and third generation L fields are 
odd. Thus, in the symmetry limit, Z2'^ forbids the dimension- five operators 
containing L2 but allows those involving Li. The same argument forbids 
operators involving L3 in the 2^^'^ symmetry limit. 

The predicted nucleon decay modes are obtained from Eq. (^Tj) by 'dress- 
ing' the two-scalar-two-fermion operators with wino exchange.^] Below, we 

^In principle, one can dress the operator udsiy by gluino exchange. Note that this op- 
erator is an SU(3) flavor singlet because of the total anti-symmetry in the color indices. 
The gluino dressing can possibly lead to the four-fermion operator (udj^sv), {ds)(ui>) or 
{su){di>). However, all these operators have mixed symmetry under flavor SU(3). Because 
of the high degeneracy among the squarks in our model, flavor SU(3) is a very good sym- 
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follow the notation of Hisano, Murayama and Yanagida |22|. The first oper- 



ator chhhs(Q2Q2)iQiLi) /M^ gives us the following four-fermion operators" 
^ = [^Xidu)isu,) - \isu)idu,)] if ice) + ficd)). (32) 



where / is the "triangle diagram factor" [P3|, a function of the wino and 
scalar masses: 



2 2 
iTLi — ml 



f(j j\ = Hi ( LI! In ^ -J 




-J V'-'l --^2 --^-Z 

and where parentheses now indicate the contraction of spinor indices. In the 
limit where M2 <^ m, this function is well approximated by / ~ M2/m^. The 
second term in (pl|), —chhhsXiQiQi)iQ2Li)/M^, gives us the four-fermion 
operators 

C = [\idu)isu,) + \isu)iue)] if ice) + fiu,d)). (34) 



All the fields in ( ^2]) and (|3^) are in the mass eigenstate basis, and terms of 
higher orders in A have been neglected. The sum of these results is given by 

= [-\isu)idu,) + \isu)iue)] ific e) + /(c, d)). (35) 

Here we have used the fact fic,d) = fiu,d) to good accuracy. Note that 
there is a precise cancellation between idu)isvg) operators in (^) and (|3^) . 
This has an important effect on the relative branching ratio between the 
charged lepton and neutrino mode, as we will discuss below. 

The ratios of decay widths can be estimated using the chiral Lagrangian 
technique |2|, |5|, |2[. We find 



Tip K^Ue) ■■ Tip K^e^) : r(n ^ K^v^) 

metry even with the squarks. Therefore, the four-fermion operators are suppressed by the 
small non-degeneracy among the squarks, and hence are negligible. There is a possibility 
that gluino dressing of operators involving third-generation fields may be important. But 
their coefficients are already as small as the one we discuss here, and they need to pick up 
smaller mixing angles and hence are negligible. 

"in the following discussion, we assume that the Cabibbo mixing originates from the 
down sector, i.e. a = 1, a' = in eq. (||) and (^. However we checked that all the results 
remain the same even when the Cabibbo mixing comes from both the down and the up 
sectors, or even solely from the up sector. 
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2mp ^ 


2 




3m B 





TUB 



F) 



3m 



iD-3F) 



B 



0.4 : 1 : 2.7. 

(36) 



where we have taken ms = 1150 MeV~ ~ m\, D = 0.81 and F = 0.44. 
We have stressed earher that the dominance of the electron over the muon 
mode is a remarkable feature of this model, which can never happen in SUSY- 
GUTs ^ In addition, what is remarkable about the result in (^) 
is that the proton's charged lepton decay mode dominates over the neutrino 
mode. This is a consequence of the cancellation of the {du){sh'e) operator 
in (P^l). The dominance of p — > f^^e^ over p — > K^Ve is rarely the case in 
SUSY-GUTs. 

Finally, we come to the overall rate. We find 



r{n ^ K^u,) 
= 4 X 10^^ 



years 



1 0.003 GeV^ 0.81 
c 



TeV" 



(37) 



e As (l + tan2/3) /(c, e) + /(c, d) 

This result includes the effect of running the dimension-five operator between 
the Planck scale and mz, (a factor of As = 0.81 in the amplitude if mt = 
175 GeV, tan/3 = 2, as{mz) = 0.12) and between mz and m„ (a factor of 0.22 
in the amplitude). In the expression above, ^ is the hadronic matrix element 
of the four-fermion operator evaluated between nucleon and kaon states; its 
exact value is rather uncertain, but is estimated to be within the range 
^ = 0.003-0.03 GeV^. If we take that M2 ~ 100 GeV and m^ ~ 700 GeV, 
and mj ~ 300 GeV, then the triangle functions /(c, e)+/(c, d) ~ (1.8 TeV)^^. 
Thus, if c = 1 and ^ = 0.003 GeV^, we obtain a mean lifetime 12.7 x lO'^^ 
years, which can be compared to the experimental bound, r(n K^Ve) > 
8.6 X 10^^ years. It is interesting to note that the coefficient 4 x 10^^ in 



would be the same in the minimal SU(5) SUSY-GUT with an extremely 
large color-triplet Higgs mass M^f^ = 10^^ GeV. Thus, the rate in our model 
is roughly comparable. Overall, the (5*3)^ symmetry gives us just enough 
suppression of dimension-five operators to evade the current bounds, so the 
model is phenomenologically viable. Since the SuperKamiokande experiment 
is expected to extend Kamiokande's current reach by another factor or 30, 
there is a very good chance that the n K^Ve mode may be seen. It is an 
exciting prediction of this model that the p K^e^ and K^v^ modes are 
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likely to be seen at the same time because their rates are close to each other, 
as we saw in eq. (|36D. 



7 Conclusions 

We have shown in this paper how to incorporate the leptons in the (5*3)^ 
model presented in Ref. [|l^. By assuming that all possible nonrenormaliz- 



able operators are generated at the Planck scale, we showed that the trans- 
formation properties of both the leptons and the flavor symmetry breaking 
fields could be uniquely determined. We then demonstrated that the phe- 
nomenological constraints from lepton flavor violation, CP violation, and 
proton decay are indeed satisfied in our model. The most striking prediction 
that emerged from our analysis is the dominance of proton decay to final 
states involving first generation lepton fields, unlike the case in SUSY GUTs. 
We showed that the ratios of decay widths for the largest modes n K^u^, 
p K^Vf., and p — > K^e^ are approximately 0.4 :: 1 :: 2.7. Given our 
estimate of the total rate, we pointed out that all three modes may be within 
the reach of the SuperKamiokande experiment and could well be discovered 
simultaneously. 
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